Numerical Evidence for a p x — ip y Paired Fractional Quantum Hall State at v — 12/5 
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We provide numerical evidence supporting a Bonderson-Slingerland (BS) non-Abelian hierarchy state as a 
candidate for the observed v = 12/5 quantum Hall plateau. We confirm the existence of a gapped incompress- 
ible v = 12/5 quantum Hall state with shift S — 2 matching that of the BS state. The exact ground-state 
of the Coulomb interaction state on the sphere is shown to have large overlap with the BS ground-state trial 
wavefunction. The analysis of the BS states is extended to hierarchical descendants of general paired states in 
the weak-pairing phase at v — 5/2. 
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Fractional quantum Hall (FQH) physics in the lowest Lan- 
dau level is well understood in terms of the Laughlin states (Hi 
and the Haldane-Halperin (HH) hierarchy states 1311, which 
can equivalently be described using Jain's composite fermion 
(CF) approach J3, H- The first appearance of an even- 
denominator fractional plateau at v = 5/2 made it clear that 
the physics of the second Landau level (2LL) could be even 
more interesting. Numerical studies J(J 0, S 13 El sup por t 
the non-Abelian p x — ip y paired Moore-Read (MR) state |1 1] 
(and its particle-hole conjugate, MR) as the correct descrip- 
tion of the v = 5/2 FQH state. At first, it seemed that this 
exceptional filling fraction was just an anomaly that appeared 
amongst other "standard" odd-denominator FQH states at v = 
7/3, 8/3, and 14/5 EH- Later, a v = 12/5 plateau also 
emerged 13, and it was numerically shown that in addition to 
the Abelian HH state, the particle-hole conjugate of the non- 
Abelian 3-clustered Read-Rezayi (RR) state il5l[l6ll is also a 
viable candidate for this filling fraction. In fact, it has been 
shown numerically that pairing/clustering is generally rele- 
vant in the 7/3 < v < 8/3 range JH. 

Recently, a non-Abelian hierarchy of states constructed 
over the v = 5/2 MR state was proposed to describe all 
the 2LL FQH states 111 811 . These Bonderson-Slingerland (BS) 
states exhibit the same pairing as the parent MR state, thus 
suggesting that the physics of the v = 5/2 "anomaly" could 
in fact be representative of all 2LL states. There has been 
much recent interest in non-Abelian FQH states due to their 
potential use for topologically protected quantum computa- 
tion 1 13 11. While the RR state can provide computationally 
universal gates from braiding alone, the BS states cannot, re- 
quiring at least one supplemental unprotected gate. Hence, 
the HH, BS, and RR v = 12/5 candidate states have vastly 
different levels of utility for quantum computation, and dis- 
covering which of these actually occur in experiments will be 
quite significant. In this Letter, we provide numerical evi- 
dence establishing the BS state as a competitive candidate at 
v = 12/5. 



The BS hierarchy states 11811 built over the MR state 
are constructed by successively condensing minimal charge 
Abelian quasiparticles and projecting them into new FQH 
states. They can be succinctly described as Ising x U(l)/f| c , 
where the coupling constant A'-matrix has A'oo = 2 cor- 
responding to the MR parent state and C is the topological 
charge spectrum. Some of these states can also be described 
using an equivalent CF type formulation 11811 . Among these 
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is the BS state with K 



, which is a candidate for 



12/5. It has the CF type ground-state wavefunction 13211 
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where Vlll is lowest Landau level projection, \n is the 
wavefunction of n filled Landau levels (n < correspond- 
ing to negative flux), xl/f^ 1 ^ is the bosonic v = 1 MR ground- 



(CF) 

state wavefunction 111 111 , and ^ 2 is the standard v = 2/3 CF 



ground-state wavefunction 
on the sphere, where 



1. This BS state has shift S = 2 



N rh = is- 1 N e -S 



(3) 



is the relation between the number of flux quanta and the 
number of electrons N e . The HH and RR states at v = 12/5, 
respectively, have 5 = 4 and —2 on the sphere. In order 
to study the validity of the BS state and to compare it with 
these other candidates for v = 12/5, we used a combination 
of powerful numerical techniques on the sphere: exact diag- 
onalization, variational Monte Carlo, and the density matrix 
renormalization group (DMRG) method. 

A necessary signature of a FQH state is the existence of a 
charge excitation gap 



A (iV ) = E N ^ +1 + En^-x - 2E N ^ 



(4) 
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FIG. 1: The charge excitation gaps for a system of N e = 14 elec- 
trons, a) A scan of gaps as a function of magnetic flux. Shifts corre- 
sponding to candidate states are labeled, b) A variation of the pseu- 
dopotential SVi around the Coulomb point at fluxes corresponding 
to v = 12/5 at shifts S = 4,2, and -2. 



at the corresponding N<j, given in Eq. (0), where E^^ is the 
ground-state energy (in units of e 2 /£o, where £o = y/ hc/eB 
is the magnetic length) for the given value of N$ fluxes. A/n 
is the energy gap of a quasihole-quasielectron pair, where n 
is the number of fundamental quasiholes produced per flux. 
(n = 2 for the HH, BS, and RR states at v = 12/5.) As 
the v = 12/5 candidate states that are being considered all 
have distinct shifts, the existence of charge gaps can be used 
to help identify which states are competitive. It is, however, 
also important to remember that finite systems can run into 
the aliasing problem, where two states with different filling 
fractions share the same value of for a given N e . 

In a recent numerical study l2lll with finite layer widths, a 
v = 12/5 state with 5 = 2 was clearly identified, with charge 
gaps given for up to N e = 14. Second Landau level flux scans 
were only performed for N e 



10 and 12 in Ref. |21|], and, 
unfortunately, at these system sizes there are aliasing conflicts 
between v = 12/5 states with 5 = 4 and 5 = —2 and the v = 
5/2 MR and v = 7/3 Laughlin (L7/3) states, respectively. 
In order to overcome these aliasing difficulties, we used the 
DMRG technique of Ref. H (see also J22|]) to study system 
sizes of up to N e = 18 electrons. 

The DMRG technique belongs to the family of variational 
methods, and includes ingredients of exact diagonalization 
and numerical renormalization group. However, no a priori 
assumptions about the physics of the variational wavefunction 
are made. The algorithm relies on a truncation of the Hilbert 




FIG. 2: The pair correlation function for the v = 12/5 state with 
S = 2 at the Coulomb point for N e = 12, 14, and 16. 



space in such a way that the loss of information is minimized. 
The resulting variational wavefunction is the best approxima- 
tion to the actual ground-state in the form of a matrix product 
state. The accuracy is completely under control and it is dic- 
tated by the number of DMRG states m kept in the truncation. 
In this work, we use up to m = 4000 states for the larger sys- 
tem sizes, giving estimated errors in the energies per electron 
on the order of 5 x 10 -5 in the worst cases. 

In Fig. QJ), we show a scan of the charge gap as a func- 
tion of magnetic flux, for N e = 14 at the Coulomb point. 
We can identify different states according to their shift, la- 
beling the v = 5/2 MR state, the v = 12/5 HH, BS, and 
RR states, and the v = 7/3 L 7 / 3 state. We find gaps for 
v = 12/5 states at S = 2 (as referenced earlier in 11811 ) and 
S = —2, but not at S = 4. In Fig. QJi), we show the behav- 
ior of the charge gap as a function of the pseudopotential Vi 
varied around the Coulomb point, for N e = 14. This exhibits 
the generally observed behavior where increasing V\ destroys 
the non-Abelian clustered states (BS and RR) and stabilizes 
the Abelian state (HH). We note that the S = 2 and -2 states 
both show a strong gap in the same range of 6V1, including at 
the Coulomb point (SVi = 0). These three v = 12/5 states 
satisfy the L 2 = condition of FQH ground-states when their 
gaps are positive, until 8V\ < —0.02. 

To further investigate the characteristics of the v = 12/5 
state with S = 2, we calculate the pair correlation func- 
tion g (r) obtained from exact diagonalization. The results 
at the Coulomb point are displayed in Fig. [2] These exhibit 
strongly damped long-distance oscillations indicative of an in- 
compressible state, further corroborating that this is indeed 
a good FQH state. We also see a slight "shoulder" at small 
r, which becomes more pronounced as SVi decreases. Such 
shoulders are present for the MR and RR states Il5l 12311 . and 
are considered a characteristic of non-Abelian clustered states. 

The preceding discussion of the spectral properties of the 
Coulomb Hamiltonian in the 2LL reveals clear evidence of the 
existence of an incompressible state at v = 12/5 with 5 = 2, 
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wavefunction as follows: 
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FIG. 3: Squared overlaps for N e — 8, 10, 12, and 14 between the 
exact diagonalization ground-state and: a) the BS ground-state wave- 
function of Eq. l[2j, and b) the BS ground-state with optimized pair- 
wavefunction of Eq. l[5}. Error bars represent the statistical uncer- 
tainty of the Monte-Carlo sampling of the overlap integral. 



consistent with the BS state. However, finding such a state 
at the same filling fraction and shift as a proposed candidate 
is still only circumstantial evidence, and more direct evidence 
is necessary to establish the BS state as an accurate descrip- 
tion. We therefore consider the overlap of the ground-state 
wavefunction of Eq. (|2| with the exact ground-state wavefunc- 
tion obtained for the same shift on the sphere. As shown in 
Fig-EtX these overlaps reach as high as 0.989(2) for N e — 8 
and remain as large as 0.83(2) for the largest system consid- 
ered (N e — 14) at SV\ = 0. Again, we manipulated the first 
pseudopotential coefficient V\ around the Coulomb potential 
of a thin 2DEG in 2LL to obtain a simple parametrization of 
the relevant interactions. The largest values of the overlap are 
obtained at slightly positive values of 5V\ ~ 0.01. The nu- 
merical evaluation of the overlap integral was undertaken by 
a Monte-Carlo sampling of O = J d{z\, . . . , zn)^*^^^ 
in position space. The composite fermion part F) of the 

3 

trial state in Eq. (|2]i was generated as a Slater determinant of 
individually projected CF orbitals 12411 at negative effective 
flux Il25ll . The rate-limiting step is the evaluation of the exact 
wavefunction, which requires calculating a number of Slater 
determinants equal to the dimension of the Hilbert-subspace 
Dl,=o projected onto fixed L z = 0. For our largest system, 
N e = 14, we have D L ^ ~ 1.9 x 10 7 . 

The MR state may be regarded as one representative of an 
entire family of weakly paired CF states @|2^]. Similarly, this 
holds true for the BS states that are derived from the MR state 
by condensation of quasiparticles. Other representatives in 
either class of paired states can be obtained explicitly by vary- 
ing the pair wavefunction In this variational approach, 
we introduce a number of parameters gk to replace the pair- 



J^.9k0k (Zi) 0- k (zj) 



L k 
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where </>k (2^) = Jf 1 V LLL {4>\ i {zj) Jj] are the projected CF or- 
bitals, and Ji = Yik^i ( Zi ~ Zk )- To obtain a paired state at 
the shift of the BS state, an effective field with —1 flux quanta 
is required for these CFs. The pair-wavefunction then acquires 
a phase of — 27r when two CFs are braided, which we denote 
as p x — ip y , or negative p-wave pairing. Our simulations are 
undertaken on the sphere, where the expansion in Eq. © in- 
volves monopole harmonics (for details, see l27ll . App. A). 
As for the paired state at v = 5/2, the number of relevant pa- 
rameters, gk, on the sphere is small [9fl - only up to 5 for the 
system sizes considered. 

Fig.Et>) shows results for the overlaps of BS states with pair 
wavefunctions optimized such as to yield maximum overlap 
with the exact ground-state at v = 12/5 and shift S = 2. 
Comparing to the results in Fig. [3^), we find that the re- 
gion of large overlaps with the exact ground-state becomes 
wider, while the overlap peaks increase significantly and shift 
to slightly higher 8V\. The overlaps now reach as high as 
0.990(2) for N e = 8 and climb to 0.92(3) for our largest sys- 
tem (JV e = 14) at 5Vi = 0.02. 

For v = 5/2, the weakly paired states are continuously 
connected |9] to a CF Fermi-liquid state (similar to the one 
at v = 1/2) at large 8V\, where the CF formulation becomes 
virtually exact. At v = 12/5, the HH state occurs at a different 
shift, so it comes as no surprise that the overlap of the BS state 
decreases for large 8V\ . We find a discontinuous drop to zero 
of the overlap in some cases (N e = 10 and 14), indicating 
level crossings in the ground-state. 

The large overlaps between the BS ground-state and the 
exact ground-state at v = 12/5 with 5 = 2, together with 
the evidence for a gapped, incompressible non-Abelian FQH 
state at this filling factor and shift, clearly establish the BS 
state as a strong candidate for the observed v = 12/5 FQH 
state 11411 . joining the ranks of HH and RR as the primary con- 
tenders. Naturally, we would like to pit these states against 
each other to see which emerges victorious. However, this 
is not so easily accomplished with numerics. For example, 
since these states have different shifts on the sphere, one can- 
not directly compare energetics, e.g. it would not be valid to 
claim the larger gap exhibited in Fig. QJ>) favors RR over BS. 
In order to compare the energetics in a somewhat meaning- 
ful way, we attempt a finite size scaling to the thermodynamic 
limit, where the shift becomes irrelevant. When comparing 
states at different shifts, we use the rescaled magnetic length 

= ^Nz/vN^iv and units of energy e 2 /£' , which compen- 



sates for finite size effects in spherical systems 112811 . 

In Fig. |U we plot the rescaled ground-state energies per 
electron at the Coulomb point for the shifts corresponding 
to the candidate v = 12/5 states, and use a least-squares 
fit to linearly extrapolate to the thermodynamic limit. Even 
though the RR state has lower energy in finite systems, the 
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FIG. 4: Finite-size scaling of ground-state energies per electron as a 
function of 1/N e for v = 12/5 states with S = 4, 2, and -2 at the 
Coulomb point. 



ground-state energies per electron in the thermodynamic limit, 
E/N e = -0.3416(5), -0.342(3), and -0.3421(5) for S = 
4, 2, and —2, respectively, are very close and within extrap- 
olation errors of each other ll33ll . so there does not appear to 
be a clear preference between them. Hence, it is likely that 
the physical effects of finite layer thickness and Landau level 
mixing will play an important role in determining which state 
is actually favored and experimentally realized, and a more 
thorough analysis of such factors is certainly warranted. 

Another way to more directly compare different states of 
the same filling fraction is to examine them on the torus, 
where all states trivially have zero shift. The numerical work 
m Ref. 1 16] examined a particle-hole symmetric system at 
v = 13/5 for A^ e = 15 and 18. The results exhibited ground- 
state degeneracy on the torus that agrees with the HH state 
for most of the parameter space, and best agrees with the RR 
state in a small region near the Coulomb point. However, the 
gap is not so large in this region, and close inspection of the 
numerics also reveals low lying states that may be identified 
as BS ground-states j29tl . This again indicates it is likely that 
the inclusion of important physical effects will be significant 
in determining which state is actually energetically favored. 
Furthermore, no scaling analysis was carried out in Ref. iflrjll . 
so, as we have shown on the sphere, it is still unclear which 
will be favored in the thermodynamic limit. 

It will be very interesting to see which state experiments 
support as the correct description of the v = 12/5 FQH 
plateau (or the so far unobserved v = 13/5 plateau). In- 
deed, it may even turn out that more than one of these states 
can be experimentally obtained by realizing different physi- 
cal regimes. Experiments that measure the electric charge of 
the fundamental quasihole will not distinguish between HH, 
BS, and RR, since these all have e/5 charged fundamental 
quasiholes. Experiments that probe scaling behavior or ther- 
mal conductance may potentially be able to distinguish be- 
tween these states 13011 . but will likely be obfuscated by non- 
universal effects. Interference experiments, however, should 



be able to unambiguously distinguish between these possibil- 
ities |31]. 
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